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— , ■ Abstract 

| In this paper we employ a "direct method" in order to obtain rank-k solu- 

£> ■ tions of any hyperbolic system of first order quasilinear differential equations 

in many dimensions. We discuss in detail the necessary and sufficient con- 
ditions for existence of these type of solutions written in terms of Riemann 
t-h ■ invariants. The most important characteristic of this approach is the introduc- 

tion of specific first order side conditions consistent with the original system of 
PDEs, leading to a generalization of the Riemann invariant method of solving 
multi-dimensional systems of PDEs. We have demonstrated the usefulness of 
our approach through several examples of hydrodynamic type systems; new 

^ (— I classes of solutions have been obtained in a closed form, 

-i— > ■ 

£h ! Resume 

S; 

J> \ Dans cet article, nous employons une "methode directe" pour obtenir des so- 

I^j ' lutions de rang k pour tout systeme hyperbolique d'equations differentielles 

quasilineaires du premier ordre en plusieurs dimensions. Nous discuterons en 
detail les conditions necessaires et sufiisantes pour l'existence de ces types de 
solutions ecrites en termes d'invariants de Riemann. La caracteristique la plus 
importante de cette approche est l'introduction de contraintes differentielles 
de premier ordre supplementaires et compatibles avec le systeme d'EDPs orig- 
inal, ce qui conduit a une genera-lisation de la methode des invariants de Rie- 
mann pour la resolution de systemes d'EDPs en plusieurs dimensions. Nous 
avons demontre l'utilite de cette approche par plusieurs exemples de systemes 
de type hydrodynamique et des nouvelles classes de solutions sont obtenues. 
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1 Introduction 

This work has been motivated by a search for new ways of constructing multiple 
Riemann waves for nonlinear hyperbolic systems. Riemann waves and their superpo- 
sitions were first studied two centuries ago in connection with differential equations 
describing a compressible isothermal gas flow, by D. Poisson ^J] and later by B. 
Riemann [20]. Since then many different approaches to this topic have been devel- 
oped by various authors with the purpose of constructing solutions to more general 
hydrodynamic-type systems of PDEs. For a classical presentation we refer reader 
to a treatise by R. Courant and D. Hilbert [2] and for a modern approach to the 
subject, see e.g. and references therein. A review of most recent develop- 

ments in this area can be found in [3 ^] . 

The task of constructing multiple Riemann waves has been approached so far 
through the method of characteristics. It relies on treating Riemann invariants as 
new independent variables (which remain constant along appropriate characteristic 
curves of the basic system). This leads to the reduction of the dimensionality of 
the initial system which has to be subjected however to the additional differential 
constraints, limiting the scope of resulting solutions. 

We propose here a new (though a very natural) way of looking at solutions ex- 
pressible in terms of Riemann invariants, namely from the point of view of their 
group invariance properties. We show that this approach (initiated in jUH2]) leads 
to the larger classes of solutions, extending beyond Riemann multiple waves. 

We are looking for the rank-k solutions of first order quasilinear hyperbolic system 
of PDEs in p independent variables x l and q unknown functions u a of the form 



We denote by U and X the spaces of dependent variables u = (u 1 , . . . ,u q ) G M 9 
and independent variables p ) G W, respectively. The functions 

are assumed to be real valued functions on U and are components of the tensor 
products A^di <g> du a on X X U. Here, we denote the partial derivatives by u° = 
diU a = du a /dx l and we adopt the convention that repeated indices are summed 
unless one of them is in a bracket. For simplicity we assume that all considered 
functions and manifolds are at least twice continuously differentiable in order to 
justify our manipulations. All our considerations have a local character. For our 
purposes it suffices to search for solutions defined on a neighborhood of the origin 
x — 0. In order to solve (Jl.lj) . we look for a map / : X — > J l (X x U) annihilating 
the contact 1-forms, i.e. 



The image of / is in a submanifold of the first jet space J 1 over X given by 
for which J 1 is equipped with coordinates x l ,u a ,uf. 




Ii = l,..., I. 





(1.2) 
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This paper is organized as follows. Section 2 contains a detailed account of 
the construction of rank-1 solutions of PDEs In section 3 we discuss the 

construction of rank-k solutions, using geometric and group invariant properties 
of the system (jl.lj) . Section 4 deals with a number of examples of hydrodynamic 
type systems which illustrate the theoretical considerations. Several new classes of 
solutions in implicit and explicit form are obtained. Section 5 contains a comparison 
of our results with the generalized method of characteristics for multi-dimensional 
systems of PDEs. 



2 The rank-1 solutions 

It is well known j2] that any hyperbolic system admits rank-1 solutions 

u = /(r), r{x, u) = Xi(u) x\ (2.1) 
where / = (/ a ) are some functions of r and a wave vector is a nonzero function 

A(u) = (Ai(u),...,Ap(u)) (2.2) 
such that 

ker (A*Aj) ± 0. (2.3) 

Solution ()2.1j) is called a Riemann wave and the scalar function r(x) is the Riemann 
invariant associated with the wave vector A. 

The function / is a solution of (jl.l|) if and only if the condition 

(a*(/)a,(/)) r , r = d £ (2.4) 

holds, i.e. if and only if /' is an element of ker (A* Aj). Note that equation (|2.4j) is an 
underdetermined system of the first order ordinary differential equations (ODEs) for 
/. The image of a solution (|2.1|) is a curve in U space defined by the map / : K — > M. q 
satisfying the set of ODEs ()2.4|) . The extent to which expresion (|2.4|) constrains the 
function / depends on the dimension of ker (A l Aj). For example, if A* Aj = then 
there is no constraint on the function / at all and no integration is involved. The 
rank-1 solutions have the following common properties : 
1. The Jacobian matrix is decomposable (in matrix notation) 

\ dr du} dr ' ^ ^ 

or equivalent ly 

^(i-I^V'a, (2.6) 



dr \ du dr 
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where we have 



du = «) e R qxp , %=( E W 1 , 

dr \dr J 



du \du a J du 

This property follows directly from differentiation of (|2.1|) . The inverses (l — f^) 

or (l — §^§£) 1 are scalar functions and are defined, since dr /du = at x = 0. From 
equations (J2.5|) or (j2.6j) . it can be noted that u(x) has rank at most equal to 1. 

2. The graph of the rank-1 solution T = {x,u(x)} is (locally) invariant under the 
linearly independent vector fields 

X a = £(u)di, a=l,...,p-l (2.8) 

acting onIx(J space. Here the vectors 

e«(«) = (ei(«),...,«(«)) T (2.9) 

satisfy the orthogonality conditions 

Ai£ = 0, a=l,...,p-l (2.10) 

for a fixed wave vector A for which ()2.3|) holds. The vector fields ()2.8|) span a 
Lie vector module g over functions on U which constitutes an infinite-dimensional 
Abelian Lie algebra. The algebra g uniquely defines a module A (over the functions 
on U) of 1-forms \i(u) dx l annihilating all elements of g. A basis of A is given by 

\ = \ i (u)dx\ £\i = Q (2.11) 

for all indices a = 1, . . . ,p — 1. The set {r = \i(u)x\ u 1 , ... , u q } is the complete set 
of invariants of the vector fields (|2.8() . 

3. It should be noted that rescaling the wave vector A produces the same solution 
due to the homogeneity of the original system (jl.lj) . 

4. Due to the orthogonality conditions (|2.10j) . together with property ()2.5j) or ()2.6|) . 
any rank-1 solution is a solution of the overdetermined system of equations composed 
of system (jl.l|) and the differential constraints 

£(«K = 0, a = l,...,p-l. (2.12) 

The side equations (|2.12jl mean that the characteristics of the vector fields ()2.8|) are 
equal to zero. 

5. One can always find nontrivial solutions of (|2.4|) if (jl.lj) is an underdetermined 
system (/ < q) or if it is properly determined (I = q) and hyperbolic. Here, a weaker 
assumption can be imposed on the system (jl.lj) . Namely, it is sufficient to require 
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that eigenvalues of the matrix (A l Aj) are real functions. 

The method of construction of rank-1 solutions to (jl.lj) can be summarized as 
follows. First, we seek a wave vector A = (Ai, . . . , X p ) such that 

rank (A^Ai) < /. (2.13) 
For each such choice of A, we look for the solutions 7° of the wave relations 

(A^A,) 7 a = 0, fji = l,...,L (2.14) 
Functions f a (r) are required to satisfy the ODEs 

f a (r) = l a (f(r)). (2.15) 

Alternatively, the system of equations (J2.4|) is linear in the variables Aj. Nonzero 
solutions Aj exist if and only if 

rank(A^(/(r))/ to (r)) < p. (2.16) 

If ()2.16|) is satisfied for some function f(r) then one can easily find Aj(r) satisfying 
equations (|2.4|) . Using u — f(r) one can define Aj (it) (not uniquely in general). If 
I < p then (|2.1fij) is identically satisfied for any function f(r) and this approach does 
not require any integration. 



3 The rank-k solutions 



This section is devoted to the construction of rank-k solutions of a multi-dimensional 
system of PDEs (II. lj) . These solutions may be considered as nonlinear superposi- 
tions of rank-1 solutions. 

Suppose that we fix k linearly independent wave vectors A 1 , . . . , \ k , 1 < k < p 
with Riemann invariant functions 

r A (x,u) = \f(u)x\ A = l,...,k. (3.1) 

The equation 

u = f(r(x,u)), r(x, u) = (r 1 (x, u), . . . , r h (x, u)) (3.2) 

then defines a unique function u(x) on a neighborhood of x — 0. The Jacobian 
matrix of (j3.2j) is given by 



\ dr du J dr 



(3.3) 
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or equivalently 



*=%(i-%%r\ (3.4) 



dr \ du dr r 
where / = {f a ), f a are arbitrary functions of r = (r ) and 

du = «) e R qxp , d J- = e R^ k , 



?L = 




dr 


\dr A , 


dr , 


fdr A \ 


du 


\du a ) 



(3-5) 



A = (\f ) e R kxp , ^-=[^—) = e R kxq . 

v ' ; du \du a J du a 

We assume here that the inverse matrices appearing in expressions (J3.3J) or (|3.4j) . 
denoted by 

\ dr du J \ du dr J 

respectively, are invertible in some neighborhood of the origin x = 0. This assump- 
tion excludes the gradient catastrophe phenomenon for the function u. 

Note that the rank of the Jacobian matrix (j3.3|) or ()3.4j) is at most equal to k. 
Hence the image of the rank-k solution is a k-dimensional submanifold S which lies 
in a submanifold of J 1 . 

If the set of vectors 

Uu) = (e a (u),...,e a (u)) T , a = l,...,p-k, (3.7) 
satisfies the orthogonality conditions 

AfC = (3.8) 
for A — 1, . . . , k, a — 1, . . . ,p — k then by virtue of ()3.3|) or (|3.4j) we have 

Q a a (x, ««) =£(u)u? = 0, a=l,...,p-k, a=l,...,q. (3.9) 

Therefore rank-k solutions, given by ([3.2)1 . are obtained from the overdetermined 
system fll . lj) subjected to differential constraints (DCs) (j3.9J) 

A">)< = 0, C(«K = 0, a = l,...,p-fc. (3.10) 

Note that the conditions (|3.9j) are more general than the one required for the exis- 
tence of Riemann k-wave solutions (see expression (|5.1|) and discussion in Section 
5). 

Let us note also that there are different approaches to the overdetermined sys- 
tem ([3.10)1 employed in different versions of Riemann invariant method for multi- 
dimensional PDEs. The essence of our approach lies in treating the problem from 
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the point of view of the conditional symmetry method (for description see e.g. [To]). 
Below we proceed with the adaptation of this method for our purpose. 

The graph of the rank-k solution T = {x, u(x)} of f)3. lUj) is invariant under the 
vector fields 

X a = £(u)di, a=l,...,p-k (3.11) 

acting on X x U C 1R P x IR 9 . The functions {r 1 , . . . , r h , u 1 , . . . , u q } constitute a 
complete set of invariants of the Abelian Lie algebra A generated by the vector 
fields fnjj. 

In order to solve the overdetermined system ()3.1Uj) we subject it to several trans- 
formations, based on the set of invariants of A, which simplify its structure consider- 
ably. To achieve this simplification we choose an appropriate system of coordinates 
onIx[/ space which allows us to rectify the vector fields X a , given by (j3.11j) . 
Next, we show how to find the invariance conditions in this system of coordinates 
which guarantee the existence of rank-k solutions in the form (|3.2j) . 

Let us assume that the k by k matrix 

n=(Af), 1< A,i<k <p (3.12) 

built from the components of the wave vectors X A is invertible. Then the linearly 
independent vector fields 

k 

x k+ i = d k+ i - 22 a^ +1 9j, 

A,j=l 

: (3.13) 

k 

x p = d p -J2 (n- 1 )!^, 

A,j=l 

have the required form (|3.11|) for which the orthogonality conditions ()3.8|) are satis- 
fied. The change of independent and dependent variables 

x 1 = r \x, u), . . . , x k = r k (x, u), x k+1 = x k+ \ . . . , x p = x p , u\ . . . , u q = u q (3.14) 

permits us to rectify the vector fields X a and get 

Xk+i = d^k+i, . . . ,X P = d SP . (3.15) 

Note that a p-dimensional submanifold is transverse to the projection (x, u) —>■ x 
at x = if and only if it is transverse to the projection (x,u) — > x at x = 0. The 
transverse p-dimensional submanifolds invariant under X^+i, ■ ■ ■ ,X P are defined by 
the implicit equation of the form 



u = f{x\...,x k ). 



(3.16) 
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Hence, expression ()3.16|) is the general integral of the invariance conditions 

u S k+i = 0, . . . , u SP = 0. (3-17) 

The system (II . lj) is subjected to the invariance conditions (|3.17j) and, when written 
in terms of new coordinates (x, u) G X x U, takes the form 

AM = °' =0 '---'^ p = ' ( 3 ' 18 ) 

or 

A^($ 2 ) _1 A = 0, ,u sk+1 =0,...,u SP = 0, (3.19) 
where the matrices $ 1 and $ 2 are given by 

« (* 2 )f=^-|^«?- < 3 - 2 °) 

The above considerations characterize geometrically the solutions of the overdeter- 
mined system f)3.10j) in the form (|3.2|) . Let us illustrate these considerations with 
some examples. 

Example 1. Let us assume that there exist k independent relations of depen- 
dence for the matrices A , . . . , A p such that the conditions 

A^Af = 0, A = l,...,k (3.21) 

hold. Suppose also that the original system (jl.l|) has the evolutionary form and 
each of the q by q matrices A 1 , . . . , A n is scalar, i.e. 

A° = I, Af = a i (u)5%, i = l,...,n (3.22) 

for some functions a 1 , . . . , a n defined on U, where p = n + 1 and for convenience we 
denote the independent variables by x = (t = n ) G X. Then the system 

(jl.l|) is particularly simple and becomes 

u t + a 1 ^)^! + . . . + a n {u)u n = 0. (3.23) 

The corresponding wave vectors 

A 1 = (-a 1 ( M ),l,0,...,0), 

: (3.24) 
A n = (-o»(u),0,...,0,l) 

are linearly independent and satisfy conditions (|3.2ip . 
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A vector function u(x, t) is a solution of (|3.23|) if and only if the vector field 
X = dt + a'^di 

defined on R n+I3+1 is tangent to the (n + l)-dimensional submanifold S = {u = 
u(x,t)} C IR n+,?+1 . The solution is thus identified with the (n + l)-dimensional 
submanifold S C IR n+,?+1 which is transverse to M n+,?+1 — > M n+1 : [x, t, u) — > (x, t) 
and is invariant under the vector field X. The functions {r(x, t, u) = (r 1 = x 1 — 
a 1 (ti)t, . . . ,r n = x n — a n (u)t), u 1 , . . . , u g } are invariants of X, such that dr 1 A ... A 
dr n A du 1 A ... A du q ^ 0. If we define i = t,u = u, then (r, t, u) are coordinates on 
]gm+<j+i an( j vector field X can be rectified 

X = d h 
The general solution is 
e = {F(r, u) = 0} 

where F : R n+I3 — > IR 9 satisfies the condition 

det /<9F<9r + 9F\ _^ Q 
\ 9r 9m (9m / 

but is otherwise arbitrary. Note that it may be assumed that 

dr 

— (x ,t ,u ) = 0, 
in which case the transversality condition is 

(dF \ 
— (x ,t ,u ) j ^ 0. 

Hence the general solution of ()3.23|) near (xq, to,Uo) is 

& = {u = f(r)}, 
where / : M. n — > M 9 is arbitrary. Thus the equation 

u = fix 1 - a l {u)t, . . . , x n - a n (u)t), (3.25) 

defines a unique function u(x,t) on a neighborhood of the point (xo,t ,u ) for any 
/. Note that 

t = 0, M (x,0) = /(x 1 ,...,x n ), 

so the function / is simply the Cauchy data on {£ = 0}. 
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Example 2. Another interesting case to consider is when the matrix Q l (or $ 2 ) 
is a scalar matrix. Then system ()3.18|) is equivalent to the quasilinear system in k 
independent variables x 1 , . . . , x k and q dependent variables u 1 , . . . ,u q . So, we have 



B A (u)u c , 



0. 



where 



B J 



A 1 A 



i \A 



If k > 2 then $ x is a scalar if and only if 

dr 1 dr k 
du ' du 







(3.26) 
(3.27) 

(3.28) 



and consequently, if and only if the vector fields A 1 , . . . , X k are constant wave vectors. 

Finally, a more general situation occurs when the matrix $ x (or $ 2 ) satisfies the 
conditions 



dx k+] 



0. 



dx v 



0. 



Then the system (13.18)) is independent of variables x k+1 , 
(EH) hold if and only if 



x 



d 2 r 



0,... 



d 2 r 



dudx k+1 ' dudx p 

Using (jSU) and we get 

du 



0. 



/,B=1 



(3.29) 
The conditions 

(3.30) 
(3.31) 



Equation (|3.31|) can be rewritten in the simpler form 

k 







£IE(n-')^f 



du 



0. 



l<l<k<i<p. 



(3.32) 



V B = 1 



Thus system (|3.18j) is independent of variables x k+1 , . . . , x p if the k by p — k matrix 
(Af), l<-B<A;<'i<pis equal to the matrix UC, where C is a constant k 
by (p — A;) matrix. In this case ()3.18|) is a system not necessarily quasilinear, in k 



U H 



independent variables x , . . . , x and q dependent variables u 

Let us now derive the neccesary and sufficient conditions for existence of solutions 
in the form (|3.2|) of the overdetermined system (|3.10|) . Substituting (|3.3|) or (|3.4|) 
into (jl.lj) yields 



Tr 



or ou I or 



0. 



(3.33) 
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or equivalently 

-l 



Tr 



or V ou or 



(3.34) 



respectively, where 

A"=(A^)eR w , /i=l,...,L (3.35) 

Given the system of PDEs (jl-ljl (i.e. functions A^u)) it follows that equations 
(|3.33|) (or (|3.34j) ) are conditions on the functions f a {r) and Xf (u) (or £*(«))• Since 
dr/du depends explicitly on x it may happen that these conditions have only trivial 
solutions (i.e. /=const) for some values of k. We discuss a set of conditions following 
from (|3.33|) or ()3.34|) which allow the system f)3. lOj) to possess the nontrivial rank-k 
solutions. 

Let g be a (p-k)-dimensional Lie vector module over C°°(X x U) with generators 
X a given by (j3.11j) . Let A be a k-dimensional module generated by k < p linearly 
independent 1-forms 

X A = X A {u)dx\ A = l,...,k 

which are annihilated by X a Eg. It is assumed here that the vector fields X a and 
X A are related by the orthogonality conditions ()3.8|) and form a basis of g and A, 
respectively. For k > 1, it is always possible to choose a basis X A of the module A 
of the form 

X A = dx iA + X? a dx ia , A = l,...,k (3.36) 

where (iA^a) is a permutation of (1, . . . ,p). Here we split the coordinates x % into 
x %A and x %a . Then from (J3.1|) we obtain the relation 

dr A dX A ■ 

°L_ = t p±x*. (3.37) 
du a du a 

Substituting ()3.37|) into equations ()3.33|) or ()3.34|) yields, respectively 

Tr^I-Q^)- 1 ^ =0, (3.38) 



or 



Tr (I -KaX^X^J =0, (3.39) 
where we use the following notation 

Q a = ^r la eW*\ K a = Va ^feR kxk , (3.40) 
or or 

f BX A \ 

Va= [d^J GRfc><9 ' < «= 1 '---.P- 1 - ( 3 - 41 ) 
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The functions r A and x %a are all independent in the neighborhood of the origin x = 0. 
The functions A M , |£, A, Q a and K a depend on r only. For these specific functions, 
equations (|3.38j) (or (|3.39p ) must be satisfied for all values of coordinates x %a . In 
order to find appropriate conditions for f{r) and X(u) let us notice that, according 
to the Cayley-Hamilton theorem, for any n by n invertible matrice M, (M _1 det M) 
is a polynomial in M of order (n — 1). Hence, one can replace equation (J3.38|) by 



Tr^Q^AJ =0, (3.42) 

where we introduce the following notation 

Q=(I-Q a; < »)- 1 det {l-Q a x ia ). 
Taking equation (J3.42)) and all its x la derivatives (with r=const) at x %a = 0, yields 

Tr\A»^-\J = 0, (3.43) 

Tr(A»Q {ai ...Q as) ^-\ S J =0, (3.44) 

where s = 1, . . . , q — 1 and (ai, . . . , a s ) denotes symmetrization over all indices in 

the bracket. A similar procedure for equation ()3.39|) yields (I3.43J) and the trace 
condition 



Tr(V^ (ai ,...,^ aa) A) =0, 



(3.45) 



where now s — 1, . . . , k — 1. 

Equation ()3.43|) represents an initial value condition on a surface in X space given 
by x %a = 0. Equations (J3.44j) (or (J3.45)0 correspond to the preservation of (J3.43|) by 
flows represented by the vector fields (J3.ll)) . Note that X a can be put into the form 



X a = d ia -\? a d A , C-Af = 0, A=l,...,k. (3.46) 
By virtue of (EOHjl . (jS3U), equations (I3~l4l or (l3~43|) take the unified form 

Tr ( A "f^f A ) =0 ' (3 ' 47) 

where either max s = q — 1 or max s = k — 1. 

The vector fields X a and the Lie module g spanned by the vector fields Xi, . . . , X p _ k 
are called the conditional symmetries and the conditional symmetry module of (jl.lj) , 
respectively if X a are Lie point symmetries of the original system (jl.lj) supplemented 
by the DCs (jSU [EJ. 
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Let us now associate the system (jl.lj) and the conditions ()3.9|) with the subvari- 
eties of the solution spaces 

B A = {(*,u«):A'2(t*)< = 0, fi = l,...,l}, 

and 

Bq = {(x,u {1) ) : COX = 0, a= a = l,...,q}, 

respectively. We have the following. 

Proposition 1. A nondegenerate first order hyperbolic system of PDEs / ti.i)) admits 
a (p-k)- dimensional Lie vector module g of conditional symmetries if and only if (p- 
k) linearly independent vector fields X\, . . . ,X p _k satisfy the conditions and 
\S.4 r /\l on some neighborhood of (x , Uq) of B — £>a H Bq. 

Proof. The vector fields X a constitute the conditional symmetry module g for the 
system (jl.lj) if they are Lie point symmetries of the overdetermined ststem (|3.1U|) . 
This means that the first prolongation of X a has to be tangent to the system (J3.1U|) . 
Hence g is a conditional symmetry module of (jl.lj) if and only if the equations 

prWjf a (A^(«) < = 0, pr^ 1 ^. =0, a = 1, . . . ,p - k (3.48) 

are satisfied on J 1 whenever the equations ()3.1()j) hold. Now we show that if the 
conditions ()3.43|) and ()3.47|) are satisfied then the symmetry criterion ()3.48|) is iden- 
tically equal to zero. 

In fact, applying the first prolongation of the vector fields X a 

pr«X a = X a + C n ^,V ° 



1 1 du<* 



to the original system (jl.lj) yields 

pr«X a (A"X) = A^i u0 ufu<* = 0, (3.49) 

whenever equations (j3.10j) hold. On the other hand, carrying out the differentiations 
of (j3.8j) gives 

A ? = (3-50) 
Comparing (j3.49p and (j3.5Up leads to 

Wi£iZ A (\f) = 0, (3.51) 
where we introduce the following notation 

= A^Z%\? . (3.52) 
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Here the new vector fields Zb are defined on U 

Z A = Z a A -^ e T U U. (3.53) 

It is convenient to write equation (|3.5H) in the equivalent form 

1i(^Z9 a ZX) = 0,/i= 1,...,/ (3.54) 

where the following notation has been used 

Oa = A^C (3.55) 

The assumption that system (jl.lj) is hyperbolic implies that there exist the real- 
valued vector fields \ A and 7^4 defined on U for which the wave relation 

(A^Af ) 7 ( A) = 0, A = l,...,k (3.56) 
is satisfied and that the U space is spanned by the linearly independent vector fields 

!A = l a A d ua eT u U. (3.57) 

Hence, one can represent the vector fields Z A through the basis generated by the 
vector fields {71, . . . , 7^}, i.e. 

Z a = h B AlB . (3.58) 

Using equations ()3.3|) and (j3.6|) we find the coefficients 

hi = m- 1 )! 

This means that the submanifold S, given by (|3.2j) . can be represented parametri- 
cally by 

a fa 

= h B Al%. (3-59) 
On the other hand, comparing ()3.3|) and ()3.58|1 gives 

< = Z a A \f. (3.60) 
Applying the invariance criterion ()3.48|) to the side conditions ()3.9|) we obtain 

w [1) x a {Qt) = c^!, .„.«;«';. (3.6i) 

The bracket [a, b] denotes antisymmetrization with respect to the indices a and b. 
By virtue of equations (j3.50j) and f)3.60|) . the right side of (|3.61|) is identically equal 
to zero. Substituting (|3.58|) into equation (|3.54|) and taking into account equation 
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()3.36|) and ()3.59|) we obtain that for any value of x G X the resulting formulae 
coincide with equations ()3.43|) and (|3.47j) . Hence, the infinitesimal symmetry cri- 
terion (j3.48j) for the overdetermined system f)3.1Uj) is identically satisfied whenever 
conditions (l3~4"3j) and (l3~47j) hold. 

The converse also holds. The assumption that the system (jl.lj) is nondegenerate 
means that it is locally solvable and takes a maximal rank at every point (xq, u^) G 
£?a- Therefore [Hj the infinitesimal symmetry criterion is a necessary and sufficient 
condition for the existence of symmetry group G of the overdetermined system 
([3. 10)1 . Since the vector fields X a form an Abelian distribution, it follows that the 
conditions ()3.43|) and ()3.47|) hold. That ends the proof since the solutions of the 
original system (jl.lj) are invariant under the Lie algebra generated by (p — k) vectors 



Note that the set of solutions of the determining equations obtained by applying 
the symmetry criterion to the overdetermined system (j3.1(Jj) is different than the 
set of solutions of the determining equations for the initial system (jl.lj) . Thus the 
system (j3.1Uj) other symmetries than the original system (jl.lj) . So, new reductions 
for the system (jl.lj) can be constructed, since each solution of system (j3.1(Jj) is a 
solution of system (jl.lj) . 

In our approach the construction of solutions of the original system (jl.lj) requires 
us to solve first the system (j3.47j) for Xf as functions of u a and then find u = f(r) by 
solving (j3.43j) . Note that the functions f*(Xf) are the functions Xf(f) pulled back 
to the surface S. The A^(/) then become functions of the parameters r 1 , . . . , r k on 
S. For simplicity of notation we denote f*(Xf) by Xf (r 1 , . . . ,r k ). 

The system composed of (j3.43j) and (j3.47j) is, in general, nonlinear. So, we cannot 
expect to solve it in a closed form, except in some particular cases. But nevertheless, 
as we show in section 4, there are physically interesting examples for which solutions 
of (j3.43j) and (j3.47j) lead to the new solutions of (jl.lj) which depend on some arbitrary 
functions. These particular solutions of (j3.43j) and (j3.47j) are obtained by expanding 
each function Xf into a polynomial in the dependent variables / and requiring that 
the coefficients of the successive powers of u a vanish. We then obtain a system of 
first order PDEs for the coefficients of the polynomials. Solving this system allows 
us to find some particular classes of solutions of the initial system (jl.lj) which can 
be constructed by applying the symmetry reduction technique. 

4 Examples of applications 

We start with considering the case of rank-2 solutions of the system (jl.lj) with 
two dependent variables (q = 2). Then (j3.47j) adopts the simplified form. 



fields Xx, ■ ■ ■ ,X p -k- 



Q.E.D. 




(4.1) 
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By virtue of ()3.43|h equation ()4.1|) can be transformed to 



Tr 



or V or V or 



A 



0. 



Using the Cayley Hamilton identity, we get the relation 
AB - I Tr AB — (B — I TrB) (A - I Tr A) 



(4.2) 



(4.3) 



for any 2 by 2 matrices A, B G 
form 



p2x2 



. Now we can rewrite ()4.2|) in the equivalent 



-Tr 
So we have 



' 0/ fdf ^df 
A M — — - ITr-^- 

or V 0r or 



(r) a - ITr i] a )\ 



df 



dot I )Tr[A M (77 a — YTTT] a )X] = 0. 



(4.4) 



(4.5) 



The rank-2 solutions require that the condition det df/dr ^ be satisfied (otherwise 
q = 2 can be reduced to q — 1). As a consequence of this, we obtain the following 
condition 



Tr[A^ a - ITr77 a )A] = 0, /i = 1, 



(4.6) 



which coincides with the result obtained earlier for this specific case |12j . One can 
look first for solutions X(u) of ()4.6|) and then find f(r) by solving ()3.43|) . Note that 
equations ()4.6|) form a system of Z(p — 2) equations for 2(p — 2) functions Xf a (u). 
This indicates that they should have solutions (say, for generic systems) if (jl.l|) is 
not overdetermined. 



Example 3. We are looking for rank-2 solutions of the (2+1) hydrodynamic 
type equations 



u\ + u j u) + A%u) 



0, i,j,k = 1,2 



(4.7) 



where A 1 are some matrix functions of u 1 and u 2 . Using the condition representing 
the tracelessness of the matrices A^w" and A^w", it is convenient to rewrite the 
system ()4.7|) in the following form 



Tr 



Tr 



1 u 1 + A 1 ^ u 2 + A\ 2 




A\ 2 



u\ u 2 
u\ u 2 



u 



u 



1 u' + A 21 



u 2 + A 22 



0. 



0. 



(4i 
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Let T be a smooth orientable surface immersed in 3-dimensional Euclidean (x, y, t) e 
X space. Suppose that T can be written in the following parametric form 

u = f(r\ r 2 ) = (^(r 1 , r 2 ), u 2 (r\ r 2 )), (4.9) 

such that the Jacobian matrix is different from zero 

Without loss of generality, it is possible to choose a basis X A of module A of the 
form 

where a A and 6 A are functions of u 1 and ii 2 to be determined. 

The rank-2 solution can be constructed from the most general solution of equa- 
tions flUJ) for X A = (-1, a A , b A ),A = 1, 2. These equat ions lead to a system of four 
PDEs with four dependent variables a A ,b A ,A = 1,2 and two independent variables 
u 1 and u 2 , 



- (A?a 2 + A 12 b 2 )^ + ((u 1 + A?y + (it 2 + A l2 )b 2 - 

+ (^a 1 + Al 2 b^ - ((u 1 + A?y + (u 2 + A^b 1 - 

- (A?a 2 + Al 2 b 2 )^ + ((^ + A?)a 2 + (it 2 + <)6 2 - 

d& 2 



(4.12) 



+ (^a 1 + Al'b 1 )^- - ((u 1 + A?)a l + (u 2 + A l2 )b x - 1) 



du 1 du 2 
(1 - (n 1 + A 21 )a 2 - (u 2 + Af)b 2 )^ + (A 2 ^a 2 + A 22 b 2 )^- 2 

- (1 - (u 1 + A 2 V - (« 2 + A ?)b 1 )^ - (Afa 1 + Afb^ = 0, 
(1 - (n 1 + Af)a 2 - (u 2 + A?)b 2 )^ + (A 21l a 2 + A 22 b 2 )^~ 2 

- (1 - (ix 1 + A 2 V - {u 2 + A 22 )b^ - (Afa 1 + Afb^ 2 = 0. 



Finally, a rank-2 solution of ()4.8|) is obtained from the explicit parametrization 
of the surface T in terms of the parameters r 1 and r 2 , by solving equations (|3.43|) 
in which X A adopt the form (J4.6j) 

((u 1 + A^a 1 + (u 2 + Al 2 )b> - + ((u 1 + A?)a 2 + (u 2 + Af )b 2 - 1)^ 

dr l • ^ - +A 2 * ) dr2 
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(4.13) 

(A 2 V + Af b 1 )^ + {Afa 2 + Afb 2 )^ + ((u 1 + A?y + {v 2 + Af)b> - 1)^ 

du 2 
dr 2 

while the quantities r 1 and r 2 are implicitly defined as functions of y, x, t by equation 

dm. 



+ ((n 1 + ^ V + (« 2 + ^ 2 )& 2 - l)£f = o, 



In the case when equation f)4.2|) does admit two linearly independent vector fields 
X A with e = —1, there exists a class of rank-2 solutions of equations (|4.12|) and 
(I4.13j) invariant under the vector fields 

X x = d t +u 1 d x , X 2 = d t + u 2 d y . (4.14) 

Following the procedure outlined in Section 3 we assume that the functions f l and 
f 2 appearing in equation ()3.2|) are linear in u 2 . Then the invariance conditions take 
the form 

x -uH = g{u l ) +u 2 h{u v ), y-u 2 t = a(u l )+u 2 b{u 1 ) } (4.15) 

where a, b, g and h are some functions of u 1 . 

One can show that if h = 0, then the solution of the system (|4.12j) . ()4.13|1 is 
defined implicitly by the relations 

x — u 1 t = g{u 1 ), y — u 2 t = a{v}) + u 2 g tU i, (4.16) 

where a and g are arbitrary functions of u 1 . Note that in this case the functions u 1 
and u 2 satisfy the following system of equations 



u\ + u\ l x + u 2 u\ + AWul - u\) + A X2 u\ = 0, 
u 2 + u x u\ + u 2 u\ + A 21 (ul - u 2 y ) + A 22 u\ = 0, 



(4.17) 



for any functions A£ of two variables u 1 and u 2 . 

If the function h of u 1 does not vanish anywhere (h ^ 0) then the rank-2 solution is 
defined implicitly by equations (|4.15J) and satisfies the following system of PDEs 



u\ + u'ul + u 2 u\ + Af[u 2 y - u\ + liu^ul + miu 1 )^] = 0, , | 

u 2 + u x u 2 x + u 2 u 2 y + Af[u 2 y -ul + l{u l )u 2 x + m^Uy] = 0, 



where A^ 2 and A 2 , 2 are any functions of two variables u 1 and u 2 . Given the functions 
I and m of u 1 , we can prescribe the functions a and b in expression (|4.15jl to find 

h = J Ibydu 1 , g = j [b — hm — ayjda 1 . (4.19) 
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For instance, consider a rank-2 solution of equations (|4.12|) and (|4.13jl invariant 
under the vector fields 

X 1 = d t + u^d X} X 2 = dt- u% (4.20) 

with the wave vectors X A which are the nonzero multiples of A 1 = (u 1 , —1,0) and 
A 2 = (w 2 , 0, —1). Then the solution is defined by the implicit relations 

x-u 1 t = g(u 1 ), 

2 i 2 ( 4 - 21 ) 

y + ut = h(u ) + u g iU i. 

and satisfies the following system of equations 

ul + v}ul + u 2 ul. + b(u x , u 2 )u], = 0, 

2 , 12, 22, / 1 2\ 1 n 

u t +u u x + u u y + c[u , u )u y — 0, 

where b and c are arbitrary functions of u l and u 2 . 

Thus, putting it all together, we see that the constructed solutions correspond 
to superpositions of two rank-1 solutions (i.e. simple waves) with local velocities 
u 1 and u 2 , respectively. According to jTUj, if we choose the initial data (t = 0) for 
the functions u 1 and u 2 sufficiently small and such that their first derivatives with 
respect to x and y will have compact and disjoint supports, then asymptotically 
there exists a finite time t = T > for which rank-2 solution decays in the exact 
way in two rank-1 solutions, being of the same type as in the initial moment. 

Example 4. Consider the overdetermined hyperbolic system in (2+1) dimensions 
(p = 3) 

du 

— — h (u ■ V)tt + ka grad a = 
dt 

da 

— + (u- V)a + k- 1 adivu = 0, (4.23) 

|^ = o, *_o, 

ox ay 

describing the nonstationary isentropic flow of a compressible ideal fluid. Here we 

/ \l/2 

use the following notations : u = (u 1 , u 2 ) is the flow velocity, a(t) = 7^ is 

the sound velocity which depends on t only, k = 2(7 — 1) _1 and 7 is the polytropic 
exponent. 
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The system (J4.23|) can be written in an equivalent form as 



Tr 



Tr 



Tr 




(4.24) 



0. 



We are interested here in the rank-2 solutions of ()4.24|) . So, we require that con- 
ditions ()3.43|) and ()3.47|) be satisfied. This demand constitutes the necessary and 
sufficient condition for the existence of a surface T written in a parametric form 
(J4.9|) for which equation (J4.10|) holds. In our case, p = q = 3 and k = 2, conditions 
(I3~l3l and (l3~T7l) become 



Trf A^A 

\ or 



1,2,3, 



(4.25) 



and 



Tr ( — ( ^ ^ \ df_^ 

\ dr \ 1 dr ^ 2 ^ 2 dr 1 J dr 



(4.26) 



respectively. Here, we assume the following basis for the wave vectors 



A; 



A X 1 A 2 
\ 2 \ 2 \ 2 

A A x A 2 



w 

2 

W 



(4.27) 



where v A and w A are some functions of u l and u 2 to be determined. The 2 by 3 
matrices r] a and the 3 by 2 matrix df /dr take the form 



dX< 



Va ~- 

df 
dr 



du c 



dXj /du 1 dXj /du 2 dXj /da 
dxfjdu 1 dxfjdu 2 dX 2 Jda 



d-i^/dr 1 du l /dr 2 
du 2 /dr x du 2 /dr 2 
da/dr 1 da/dr 2 



1,2 



(4.28) 
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Equations (|4.25|) lead to the following differential conditions 

du 1 du 1 , i , 9 .du 2 . i , n.du 2 9 . da da. 
d^ + d^ + {u kav ) d^ + {U ~ kaw ) W + <9~H~ + d^ ] = °' 

,9m 1 ,9m 1 i. i du 2 i<9u 2 . . 9 i , 9n da 
v 1 — + w 1 — + u (v ~ttt + w + (uV + kav 2 ) — 

+ (uW + ta 2 )|^ = 0, (4.29) 
dr 2 

, , 9 du 1 9<9w 1 \ / , 2 in du 2 

div^ da da 

- (a - w 2 ku 1 ) — + (av 1 + kv 2 u 2 )— + (aw 1 + kw 2 u 2 )— = 0. 

dr 2 dr 1 dr 2 

Assuming that we have found v A and w A as functions of u 1 and u 2 , we have to solve 
f!4.26|) for the unknown functions u 1 and u 2 in terms of r 1 and r 2 . The resulting 
expressions in the equations (j4.26|) are rather complicated, hence we omit them here. 
Various rank-2 solutions are determined by a specification of functions v A and w A in 
terms of u 1 and u 2 . By way of illustration we show how to obtain a solution which 
depends on one arbitrary function of two variables. 

Let us suppose that we are interested in the rank-2 solutions invariant under the 
vector fields 

X 1 = d t + u l d x , X 2 = d t + u 2 d y . (4.30) 

So, the functions r 1 = x — uH and r 2 = y — u 2 t are the Riemann invariants of 
these vector fields. Under this assumption, equations ()4.25j) and ()4.26|) can be easily 
solved to obtain the Jacobian matrix 

(43l) 

which has the characteristic polynomial with constant coefficients. This means that 
the trace and determinant of J are constant, 

(i) uli + u 2 2 = 2Ci, 

/.. 212' C a ,C 2 GR (4.32) 

{%%) M r iM r 2 — U r iU r i = C 2 . 

The trace condition (|4.32f i)) implies that there exists a function h of r 1 and r 2 such 
that the conditions 



u 



1 = C 1 r 1 + h^, u 2 = C 1 r 2 -h r i, (4.33) 



hold. The determinant condition f)4.32( ii)) requires that the function Z^r 1 ,?" 2 ) sat- 
isfies the Monge- Ampere equation 

h r i r ih r 2 r 2 - h 2 rlr2 =C, Cel. (4.34) 
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Hence, the general integral of the system ()4.23|) has the implicit form defined by the 
relations between the variables t, x, y, u l and u 2 



dh 



II. 



1 = Ci(x - uH) + —(a; - uH, y - uH) 



dr 



u 2 = d(y - uh) + ^(x - uH, y - uH), ( 4 - 35 ) 

a = a ((1 + Cxtf + Ct 2 ) ~ l/k , a eR 

where the function h obeys ([4.34)1 . 

Note that the Gaussian curvature K expressed in curvilinear coordinates (t, r 1 , r 2 ) G 
R 3 of the surface S = {t = ^(r 1 ^ 2 )} is not constant and is given by 

For example, a particular nontrivial class of solution of ()4.23|) can be obtained if we 
assume that C — 0. In this case the general solution of ()4.23|) depends on three 
parameters, ao, C\, m G K and takes the form 



u 



1 = d(x-u l t) + (1 -m) 



i ,\ rn 
X — U t x 



y — u 2 t 



u 2 = Ci(y — u 2 t) — m 
a 



2, \ i—m 

V~ ut \ ; (4.37) 
x — uH J 



a(t) 



;i + c 1 t) 2 / fc " 



Note that if C = and C\ — then the Jacobian matrix J is nilpotent and the 
divergence of the vector u is equal to zero. Then the expression 

x-uH\ 

57 , a = a , 

y — uH ) 

2 (y-uH 
u = — m — 

\x — it 1 1 

defines a solution -u = (-u^-u 2 ) to incompressible Euler equations 

Ut + (u • V)u = 0, divu = 0, a = a . (4.39) 
Note that for m = 2, the explicit form of (j4.38|) is 



-i ~y l - 2tx ± V|/ 2 + 4ta; 

W = 7T- 1 

2t 2 

2 y =F a/|/ 2 + 4tx 
w = . 



(4.40) 
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Example 5. Now let us consider a more general case of to the autonomous system 
(I4.23J) in p = n + 1 independent (t,x l ) G X and q = n + 1 dependent (a,u % ) G U 
variables. We look for the rank-k solutions, when k = n. The change of variables in 
the system (|4.23|) under the point transformation 

i=t, x 1 = x 1 - uH, . . . ,x n = x n - u n t, a = a, u = u (4.41) 

leads to the following system 
Du 



Dt °< 

Da Da , 

0, + aTr B~ —— = 0, a^O 



(4.42) 



Dx Dt 
where the total derivatives are denoted by 

D d d D d , d . . 

= TT + MfTT-, — r = — — r + u\.j — — r , 7 = 1, . . . , 71 (4.43) 

<9t *<9m 1 ' D& 3x3 x du 1 J y ' 
and the n by n nonsigular matrix S has the form 

B = I + t^. (4.44) 
ax 

The general solution of the first equation in ()4.42j) is 

u = f(x), x = (x 1 , . . . , x n ) (4.45) 

for some function / : M. n — > M. n . The second equation in (|4.42|) can be written in an 
equivalent form 

|L(ln|a(t)| fc ) +Tr[(I + iDf(x))- 1 Df(x)} = 0, (4.46) 
where the Jacobian matrix is denoted by 

Df(x) = ^(x). (4.47) 
Differentiation of equation ()4.46|) with respect to x yields 

— (lndet(I + tD/(x))) = (4.48) 
with general solution 

det (I + iDf(x)) = a(x)(3(t) (4.49) 
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for some functions a : R n — > R and (3 : R — > R. Evaluating ()4.49|) at the initial data 
t = implies a(x) = /?(0) _1 . Therefore 

det(I + Wf(x)) = ^ y (4.50) 

So, we have 
<9 

— det (I + tDf(x)) = 0. (4.51) 
ax 

Now, let us write the determinant in the form of the characteristic polynomial 

det (I + iD f (x)) = i n P n (e,x), e = i (4.52) 

where 

det (el + = e n + p n -i(z)£ n_1 + • • • + Pi(x)e + p (x). (4.53) 

Equation (|4.51|) holds if and only if the coefficients of the characteristic polynomial 
Po, ■ ■ ■ ,p n -i are constants. So, equation (|4.46j) implies that 

-- In \a{t) | k + -- In | det (I + Wf{x))\ = 0. 
Then we have, 

^ (\a(t)\ k det (I + Wf(x)))=0. (4.54) 

Solving equation ()4.54j) we obtain 

a(t) = 7 (det (I + iD f (x)))~ 1/k , 0^ 7 gM. (4.55) 
Thus, the general solution of system ()4.23|) is 

u = f{x), a{t) = 7 [1 +p n - 1 i + . . .+ Po iT 1/k , (4-56) 

for any differentiable function / : IR n — > R™ and takes the form of a constant 
characteristic polynomial on the Cauchy data t = 

P n (e) = e n + pn-xe^ 1 + ...+ Pl e + p . (4.57) 
Note that the function a is constant if and only if 

P n (e) = e n . (4.58) 
This fact holds if and only if the Jacobian matrix Df(x) is nilpotent. 
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Note that in the particular case when p = 2, the general explicit solution of ()4.23|) 
is given by 

u(x,t) = ((3 + ax)(l + at)~\ a{t) = 7(1 + at)~ 1/k , a,P,^eR. (4.59) 

An extension of this solution to the (n+l)-dimensional space X is as follows 

u{x,t) = (I + ta) _1 (/3 + ax), a(t) = 7(det (I + at))" 1//fc , (4.60) 

where (3 is any constant n-component vector and a is any n by n constant matrix. 
In this case the Jacobian matrix is constant 

Df(x) = a (4.61) 

for any iel" 

Finally, a similar computation can be performed for the case in which the vector 
function u = (m 1 ,?/ 2 ,?/ 3 ) satisfies the overdetermined system ([4.39)1 . In the above 
notation, an invariant solution under the vector fields 

X a = d t + u a d (a) , a = 1,2, 3 (4.62) 
is given by u = f{x) and is a divergence free solution 

divw = (4.63) 
if and only if the trace condition 



holds. But 

Du _ dB 

m ~ ~dt' 
Therefore div u = if and only if 



(4.65) 



TrLB-^^O, (4.66) 
or equivalently, if and only if 

J^(det£) = 0. (4.67) 

This means that the Jacobian matrix Df(x) has to be a nilpotent one and takes the 
form 



/i 2 fh 

P s -f, 
P s -I 



3 



Df{x) = I p s -p s I , (4.68) 

i'2 f2 
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where f l is an arbitrary function of two variables x 2 and x 3 and f 2 is an arbitrary 
function of one variable s = x 2 — x 3 . Note that if j\ 2 7^ f%3 then the Jacobian matrix 
Df(x) has rank 2 (otherwise f 1 is any function of s and Df(x) has rank 1). As a 
consequence, the matrix B has the form 

/I //,'■ //.'< \ 
£ = 1 + t/ s 2 -t/ 2 , det B = 1. (4.69) 

So, the condition ()4.67|) is identically satisfied. Thus, the general solution of system 
(I4.H9J) is implicitly defined by the equations 

u 1 = f 1 ( x 2 -tf 2 ( x 2 -x 3 ),x 3 -tf(x 2 -x 3 )) , u 2 = u 3 = f 2 (x 2 -x 3 ) } a = a , (4.70) 

where the functions f 1 and f 2 are arbitrary functions of their arguments. Equations 
()4.70|) define a rank-2 solution but, according to the formula for the corresponding 
principle ^Oj, it is not a double Riemann wave. 

Obviously, other choices of the wave vectors \ A (and the related vector fields 
X a ) lead to different classes of solutions. The problem of the classification of these 
solutions remains still open but some results are known (see e.g. the functorial 
properties of systems of equations determining Riemann invariants 11). 



5 Conclusions 

In this paper we have developed a new method which serves as a tool for constructing 
rank-k solutions of multi-dimensional hyperbolic systems including Riemann waves 
and their superpositions. The most significant characteristic of this approach is that 
it allows us to construct regular algorithms for finding solutions written in terms 
of Riemann invariants. Moreover, this approach does not refer to any additional 
considerations, but proceeds directly from the given system of PDEs. 

Riemann waves and their superposition described by multi-dimensional hyper- 
bolic systems have been studied so far only in the context of the generalized method 
of characteristics (GMC) The essence of this method can be summarized 

as follows. It requires the supplementation of the original system of PDEs (jl.lj) with 
additional differential constraints for which all first derivatives are decomposable in 
the following form 

r)n a k 

7rr(*) = £^(*M(i0tf(«), (5-i) 

OX A=l 

where 

(A"»Af) 7(A) = 0, A = l,...,k 
rank(A^(w)Af) < l 
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Here, £ A + are treated as arbitrary scalar functions of x and we assume that 
the vector fields {71, . . . ,7fc} are locally linearly independent. The necessary and 
sufficient conditions for the existence of k-wave solutoins (when k > 1) of the system 
(|5.1|) in terms of Riemann invariants impose some additional differential conditions 
on the wave vectors \ A and the corresponding vector fields 74, namely ^H] 

b/A, 7b] e span{7 A , 7^}, 

£ 7 ^A A G span{A A , X B }, VA ^ B = 1, . . . , k, { ' } 

where [7,4,7s] denotes the commutator of the vector fields and 75, ^ 7s denotes 
the Lie derivatives along the vector fields 7#. 

Due to the homogeneity of the wave relation ()5.2j) we can choose, without loss 
of generality, a holonomic system for the fields {71, . . . ,7^} by requiring a proper 
length for each vector 7^ such that 

[7A,7fl] = 0, VAj:B = l,...,k. (5.4) 

It determines a k-dimensional submanifold S C U obtained by solving the system 
of PDEs 

a fa 

^i = r A (f\...j k ) (5.5) 

with solution ir : F — > U defined by 

tt : (r\ . . . , r k ) - (f\r\ . . . , r k ), . . . , f*{r\ r k )) . (5.6) 

The wave vectors \ A are pulled back to the submanifold S and then X A become 
functions of the parameters r 1 ,...,r fe . Consequently, the requirements (|5.1|) and 
()5.3j) take the form 

— (x)=£ A (x)\ A (r\...,r k ), (5.7) 
d\ A 

— r G span{A A ,A s }, ^ S = 1, . . . , k (5.8) 
or" 

respectively. It has been shown ^H] that the conditions (|5.5|) and ()5.8|) ensure 
that the set of solutions of system (jl.lj) subjected to (|5.1jl . depends on A; arbitrary 
functions of one variable. It has also been proved J7j that all solutions, i.e. the 
general integral of the system (|5.7|) under conditions (|5.8|) can be obtained by solving, 
with respect to the variables r 1 , . . . , r k , the system in implicit form 

Af(r\...,rV = ^V>- ••>*•*), (5.9) 

where ip A are arbitrary functionnally independent differentiable functions of k vari- 
ables r 1 , . . . ,r k . Note that solutions of ()5.7|) are constant on (p — A;)-dimensional 
hyperplanes perpendicular to wave vectors A A satisfying conditions (|5.8|) . 
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As one can observe, both methods discussed here exploit the invariance properties 
of the initial system of equations In the GMC, they have the purely geometric 

character for which a form of solution is postulated by subjecting the original system 
to the side conditions (|5.1|) . In contrast, in the case of the approach proposed 
here we augment the system (Jl.lj) by differential constraints (|3.9|) . 

There are, however, at least two basic differences between the GMC and our 
proposed approach. Riemann multiple waves defined from (j5.1|) . ()5.5|) and ()5.8|) 
constitute a more limited class of solutions than the rank-k solutions postulated 
by our approach. This difference results from the fact that the scalar functions £ A 
appearing in expression (|5.1jl (which describe the profile of simple waves entering 
into a superposition) are substituted in our case (see expressions f!3.3|) or ()3.4|) ) with 
a q by q or k by k matrix $ x or $ 2 , respectively. This situation consequently allows 
for a much broader range of initial data. The second difference consists in fact 
that the restrictions (J5.5|) and (|5.8|) on the vector fields ^a and X A , ensuring the 
solvability of the problem by GMC, are not necessary in our approach. This makes 
it possible for us to consider more general configurations of Riemann waves entering 
into an interaction. 

A number of different attempts to generalize the Riemann invariants method and 
its various applications can be found in the recent literature on the subject (see e.g. 
jH El El El 122]) • For instance, the nonlinear k- waves superposition u = fir 1 , . . . , r k ) 
described in [TH] can be regarded as dispersionless analogues to "n-gap solution" of 
(II. 1J) which require the resolution of a set of commuting diagonal systems for the 
Riemann invariants r A , i.e. 

rii = ^)(r)r^, A = 1, . . . , k, i^j = l,...,p. (5.10) 

That specific technique involves differential constraints on the functions \i A - of the 
form [22] 

fyf'iti) ®3 U i(j) ... a / D , , /cii\ 

A A = B B » A ^ B = 1, . . . ,k, (5.11) 

Ni) 

no summation. As in the case of Riemann k-waves if the system (|5.11|) is satisfied 
for the functions \i A , then the general integral of the system ()5.10j) can be obtained 
by solving system ()5.9j) with respect to the variables r 1 , . . . , r k . 

In contrast, our proposed approach does not require the use of differential equa- 
tions (|5.10j) and therefore does not impose constraints on the functions fif^ when the 
1-forms X A are linearly independent and k < p. 

However, if one removes these assumptions and X A can be linearly dependent and 
k > p then the approach presented in T\ is a valuable one and provides an effective 
tool for classification criterion of integrable systems. 

In order to verify the efficiency of our approach we have used it for constructing 
rank-2 solutions of several examples of hydrodynamic type systems. The proposed 
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approach proved to be a useful tool in the case of multi-dimensional hydrodynamic 
type systems (jl.ljl . since it leads to new interesting solutions. 

The examples illustrating our method clearly demonstrate its usefulness as it 
has produced several new and interesting results. Let us note that the outlined 
approach to rank-k solutions lends itself to numerous potential applications which 
arise in physics, chemistry and biology. It has to be stressed that for many physical 
systems, (e.g. nonlinear field equations, Einstein's equations of general relativity 
and the equations of continuous media, etc) there have been very few, if any, known 
examples of rank-k solutions. The approach proposed here offers a new and promis- 
ing way to investigate and construct such type of solutions. 
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